Ghost free Massive Gravity in the Stiickelberg language 
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Massive Gravity in four dimensions has been shown to be free of the Boulware-Deser (BD) ghost 
in the ADM language for a specific choice of mass terms. We show here how this is consistent with 
the Stiickelberg language beyond the decoupling limit, and how the constraint required to remove 
the BD ghost arises in this framework nonperturbatively, without the use of field redefinitions. We 
emphasize a subtlety in obtaining this constraint, that has been overlooked in previous literature. In 
both the ADM and Stiickelberg formalisms the vanishing of the determinant of a Hessian guarantees 
the absence of the BD ghost. 



INTRODUCTION 

In recent years there has been considerable progress 
in the construction of massive gravity theories in 4 di- 
mensions, [1-9], culminating in the proposal of [10, 11] 
of a two parameter family of models which are ghost free 
in the decoupling limit and were explicitly shown to be 
ghost free in the full theory to fourth order, as well as 
completely nonlinear ly in specific cases. Furthermore, it 
was shown in [12] that the fourth order calculations of 
[11] could be generalized to all orders giving a full proof 
in the ADM formalism that this class of theories is ghost 
free. 

The central problem faced in any construction of mas- 
sive gravity is that a naive counting of degrees of freedom 
would lead to six modes. This can be seen in a number 
of ways, but in the ADM language it arises because the 
3-metric has 6 components. In the unitary gauge of a 
massive theory diffeomorphism invariance is lost because 
of the mass term and so the lapse and shift remain as 
non-dynamical degrees of freedom which in general do 
not enforce constraints. As such they can be integrated 
out, giving a theory of 6 propagating degrees of free- 
dom. Since a massive theory should only have 5, the 
additional 6th mode is inevitably a ghost and is known 
as the Boulware-Deser (BD) ghost [4]. 

The resolution to this problem is to construct the mass 
term in such a way that a constraint remains, which 
serves to project out nonlinearly the BD ghost. In gen- 
eral the mass term in the ADM form can be expressed 
as 



On 



A^xCra{N,N\h,j). 



(1) 



The existence of a constraint can be characterized by 
the statement of the vanishing of the determinant of the 
Hessian det['Ha6] — where the 4x4 Hessian T-Lab is 
constructed as 



T~Lab — 



_d d_ 

ON'' dNb 



{Cra{N,N\h,j)) , 



(2) 



with the shorthand N'^ — N. Whenever the determinant 
vanishes it implies that it is not possible to solve the equa- 
tions of motion to determine the lapse and shift directly 
and so at least one constraint must remain. In practice 
we require only one additional constraint from this pro- 
cess and so an equivalent way of stating this is that once 
we have solved the algebraic equations for the shift and 
substituted back into the action, the resulting action be- 
comes linear in N and so N enforces a new version of the 
Hamiltonian constraint. When this calculation is done 
with the phase space form of the action, it is necessary 
to generate two constraints in phase space to remove one 
physical degree of freedom. The second constraint simply 
comes from the commutator of the first constraint with 
the full Hamiltonian, as it does in the linearized theory. 

These arguments describe how in principle the ghost 
can be removed in a massive theory of gravity in the 
ADM language. Fortunately it has recently been shown 
in [12] that the two parameter family described in [11] 
satisfies all the above requirements, namely that once 
the shift is integrated out, the resulting Hamiltonian is 
linear in the lapse. This is consistent with the argument 
of [11] to fourth order, and the all orders argument of 
[11] in a two dimensional model. 



STUCKELBERG FIELDS 

The Stiickelberg trick has been successful in determin- 
ing the form of the mass term that will inevitably be 
ghost free in the decoupling limit. However one can see 
that away from the decoupling limit, a problem appears 
to arise in perturbation theory seemingly hinting towards 
the presence of a ghost. The result would be in clear 
contradiction with the ADM analysis, and so it behoves 
us to ascertain where the discrepancy comes from. The 
essence of the Stiickelberg trick is to write the mass term 
in a gauge invariant manner by defining the covariant 



tensor 



H^j.u = .g^j/ - riabdf_c(f>°-di,(j)'' 



(3) 



In unitary gauge (/)" — x°- and so H^^ = g^^, — rj^^, is just 
the metric perturbation. For this to transform as a tensor 
it is crucial that the 0" transform as scalars under dif- 
feomorphisms even though they are vectors under global 
Lorcntz transformations. This point makes it clear that 
the global Lorentz symmetry is quite separate from the 
local diffeomorphism invariance. Since this transforms 
as a tensor we can construct the mass term out of scalar 
contractions of H^^ = g^^Hau- The key realization of 
[11] was that the mass term is inevitably a finite order 
polynomial in the tensor /C^^ defined such that 



2IC^-IC^IC:=Hi:. 



(4) 



The allowed mass terms derived in [11] can be shown to 
arise as the expansion of the determinant ^/— 5 det [g^ + 
A/C(^] to fourth order in A [13 15]. With these terms it 
was shown that no ghosts are present in the decoupling 
limit. 



Stiickelberg trick in the decoupling limit 

With the Stiickelberg fields the absence of a ghost in 
the decoupling limit has traditionally been determined in 
one of two ways. One is that it must be possible to write 
the action in such a way that after integration by parts it 
is independent of (f>^ . If this is the case then varying with 
respect to (jp will enforce a constraint which removes the 
BD ghost. Equivalently it has become conventional to 
introduce a further Stiickelberg field n by means of the 
decomposition (p'^ — x"- + V"" + d""K. With this definition 
it is clear that tt is a scalar under global Lorentz trans- 
formations. Furthermore tt captures the physical degree 
of freedom of the helicity zero mode of the massive gravi- 
ton. An equivalent way of stating the problem is that the 
equations of motion for tt remain second order in time. 
This amounts to to requiring that up to total derivatives 
the action can be written in such way that no derivatives 
higher than tt" enter which is equivalent to the previous 
criterion. 



(Apparent) Failure of the Stiickelberg trick beyond 
the decoupling limit 

Although the Stiickelberg language has been success- 
fully used to show the absence of ghosts in the decoupling 
limit, this same logic has been applied in several papers 
[16, 17] to argue that there is a ghost in massive grav- 
ity since it can be shown that when we go away from 
the decoupling limit terms arise which are quadratic and 



higher in (jfi in the action, which would lead to the equa- 
tion obtained by varying the action with respect to (/)" 
being dynamical and thus removing this constraint. 

However it is easy to see that while the Stiickelberg 
field (p"- can be treated as a vector in the decoupling limit, 
this treatment no longer holds beyond the decoupling 
limit. (/)" is intrinsically a set of four scalars, and only 
in the decoupling limit where tt transforms as a scalar 
under diffeomorphisms is the previous counting appro- 
priate. Beyond the decoupling limit, there is no unique 
way to identify the scalar tt, and the previous trick cannot 
be used. Instead one should carefully check whether the 
system propagates four independent Stiickelberg scalar 
fields as naively anticipated (in which case the BD ghost 
would survive) or whether there exist a constraint and 
not all four Stiickelberg fields are dynamical. 

Fortunately, it is easy to see that the constraint 
found in the ADM language does indeed survive in the 
Stiickelberg formalism. There does exist an equation 
which is obtained from the equations of motion in which 
neither (p'^ nor (p^ enters. It is simply that this is not the 
equation obtained by varying the action with respect to 
(fp but is rather a linear combination of the equations 
obtained by varying with respect to 0". In other words 
if we denote by £a the equation obtained by varying the 
Lagrangian with respect to 0" 



^.^44^0 



(5) 



we find that although Sq contains terms with (p^, the 
dynamics of the four scalars are not independent. More 
precisely, upon inspection, the matrix 






(6) 



has zero determinant and is not invcrtiblc. In particular 
there exists a linear combination 



C = d'Ea = 



(7) 



for which C contains no higher than first order time 
derivatives of (jf" . Stated in different terms, the Lagrange 
multiplier for the constraint is not (p^ but is rather de- 
termined by the direction in field space set by 5(p'^ oc c". 
This means that if we perform the helicity decomposition 
(^" = x" + V"' + 9"7r, it would appear as if the equation 
of motion for tt were higher derivative. However these 
higher derivatives would be removed by using the well 
defined equations of motion for the V"' . Another way of 
stating this is that at the level of the action it is possible 
to perform field redefinitions to remove all terms with tt. 
Whilst this improves the understanding of perturbation 
theory, it is not necessary to perform these field redef- 
initions. The presence of a constraint can equivalently 
be shown nonlinearly in the original variables using the 
equations of motion. 



We shall show how this works explicitly in the next 
section, however it is worth pausing to ask why the con- 
straint does not arise in the normal manner. The reason 
can be traced to the fact that whilst (/>" is a vector under 
global Lorentz transformations, it is a set of four scalars 
under diffeomorphisms. The normal logic that varying 
with respect to (j)^ should determine the missing con- 
straint would naively apply if 0" were a vector under dif- 
feomorphisms. For instance this would be an appropriate 
logic for a massive gauge theory. In the decoupling limit 
this distinction becomes irrelevant, but in the presence of 
gravity it becomes extremely important. By working in 
unitary gauge, the ADM phase space analysis bypasses 
this peculiarity and gives a clear counting of degrees of 
freedom. Nevertheless as we shall show in the next sec- 
tion the counting, when done properly, is identical in the 
Stiickelberg picture. 



EXPLICIT EXAMPLE 

Locally Inertial Frame 

In this section we will focus solely on the mass term 
of the Lagrangian. Even though the absence of ghost 
in the ADM formalism has been proven for the full two 
parameter family of mass terms [12], we focus in what 
follows on the specific potential /C^j^/C — /C^, although 
the same results holds for any mass term in this class 
of models. Since no derivatives of the metric appear in 
this term, it is convenient to make a local transforma- 
tion so as to be in the Locally Inertial Frame (LIF) and 
set the metric locally to 5^^ = 77^^ and ChristofFels to 
zero. Including the metric only alleviates the problem 
to be described since the metric would otherwise carry 
some of the derivatives and in turn reduces the number 
of derivatives on the scalar fields. By working in the LIF, 
we therefore ensure to keep the largest possible number 
of derivatives on the Stiickelberg fields. We can then 
evaluate the equations of motion for the scalars in the 
LIF and these must be well defined. To start with, for 
pedagogical simplicity only, let us consider the specific 
situation where the four Stiickelberg scalar fields are set 
to (j)°- = (0°(i,x),(/)^(i,x),O,O) (or to be more precise, we 
focus on a 2-dimensional subclass of models). For that 
specific configuration, the mass term is then (up to total 
derivatives and an overall normalization) 



Lm = /C - K,*^,, /C 



= ^{d,r? F^^ 



(8) 
(9) 



with F^j_u — df_i(j)u — du<j}fi- It is already clear that due to 
the term df^cj)^ , the Lagrangian is not linear in (jp , which 
is potentially worrisome. To make this more explicit, 
one could split the Stiickelberg in a similar way as what 



is done in the decoupling limit, i.e. </>" = a;" + V"' and 
expand the mass in powers of V"' to obtain 



r 



1 



-F,.(l 



\d.v^ 



(10) 



-^^lu + -Ad^v^f + ' 



1 , 1 

— T + -( 

with now T^j^u ~ df^Vi, — d^Vfj_. Now it is straightforward 
to see that if V''^ was the O"^ component of a vector, the 
last term (9^1^^)^ would be problematic as it makes V'^ 
dynamical. Saying it differently, if we were to write V'^ = 
■n then the equations of motion would involve a term of 
the form tt, which would naively manifest the presence 
of a ghost. If this conclusion was true, it would suggest 
an inconsistency with the ADM analysis [11, 12] which 
showed the absence of ghosts in the theory. However both 
results agree when pushing the Stiickelberg analysis a 
step further and properly counting the number of degrees 
of freedom, as we show explicitly below. 

As emphasized earlier, (j)^ is not a vector but rather a 
combination of scalars, so beyond the decoupling limit, 
it is not straightforward to identify the helicity zero com- 
ponent. In particular just naively performing a decom- 
position of the form (/)" = x" -|- V'^ + d°"iT, tt alone does 
not capture the helicity zero mode. Instead, to check the 
presence or not of the 6**^ BD mode, one should directly 
count how many fields are physically propagating. To 
avoid the ghost, only three out of the four Stiickelberg 
fields should be dynamical while the last one should prop- 
agate a constraint. However since (j)^ is not a vector, 
there is no reason why 4>^ (or V''^) should be the non- 
dynamical degree of freedom. As it happens here, it is 
actually a combination of all the Stiickelberg fields which 
generates the constraint and reduces the number of de- 
grees of freedom. 

To see this explicitly, and in a completely non- 
perturbative way, let us return to the action for the two 
Stiickelberg fields (fp-'^ . The field equations of motion are 

£u - a^n + ^5,9^(/)^ (11) 



2 



d,{{do,rY-Fip) (I 



(5„r) 



,a^2 



f2 



-6^,d,r - n ) = . 



It is clear that the equation for cffi involves a dangerous 
term with (fp and that degree of freedoms seems a priori 
dynamical. However upon look, the kinetic matrix Aab 
has a vanishing eigenvalue, 



'^ab *^r, 



(/)0' + </.! d^q}- 



(d^r) 



.(12) 



So even though all the equations of motion seem dynami- 
cal, they arc not independent. This can be made manifest 
by constructing the following linear combination 

c - a^(/)^£o-n% (13) 



which involves no double time derivative term. There is 
therefore a constraint to be satisfied which prevents all 
the Stiickelberg fields to be independent propagating de- 
grees of freedom. This constraint is precisely analogous 
to that obtained in the ADM language in 2 dimensions, 
[11] and then recently completely generally in four di- 
mensions [12]. 

Saying it differently, if one were to solve the equation 
for (f)^ using £i and plug that solution back into £q, we 
would then not be able to solve for 0° but would instead 
obtain the constraint mentioned above. 

Or yet in other words, if we were to write V'^ = tt, this 
would correspond to a perfectly well defined equation of 
motion for tt. 



In the rest of the paper we focus on the four- 
dimensional scenario, where the counting goes as fol- 
lows, the massless graviton carries two degrees of free- 
dom, which added to the four Stiickelberg fields would 
lead to the usual five helicity modes of the graviton plus 
a sixth BD ghost. As we will see below, the theory 
of massive gravity presented here has a primary con- 
straint which then generates a secondary constraint both 
of which ensure the absence of the BD ghost. In two di- 
mensions there is are accidental tertiary and quaternary 
constraints which are necessary to remove all the degrees 
of freedom of the Stiickelberg fields. 



GENERAL FORMALISM 



Including Gravity 

It is instructive to follow this 2d toy-model further 
and check the algebra of the whole constrained sys- 
tem. For this we include gravity ds^ =- —N'^dt^ + 
7ii (dx^ + N^dt) . The Lagrangian is then of the form 



Cm = 2A^V7rr -i + V{'D-^-)iv+^+) 



-i + l.(v^4>-){v+^+) + ^ 



(14) 
(15) 



where A is introduced as a Lagrange multiplier and we 
use light-cone coordinates, 



6± - 0O±0i 
1 



(16) 



1 



V± = —= d^±- [dt " N^d^] . (17) 

The associate conjugate momenta are then 

P± = ±^ (P^0T) ^ (18) 

giving the following expression for the Hamiltonian 

n = — ^\ - 2-111 + p+d^(t)+ -p-d^,^-] (19) 




where A = AA^. From this expression it is clear that N , 
N^ and A are three Lagrange multipliers for three con- 
straints, and that the Hamiltonian is pure constraint. 
Consistency of the constraint algebra will generate a 
fourth constraint, which is sufficient to remove all phys- 
ical degrees of freedom. This is consistent with the ex- 
pectation that a massive graviton in two dimensions has 
no degrees of freedom. The existence of additional con- 
straints other than the two expected is an accident of two 
dimensions and does not follow in higher dimensions. 



The previous configuration shows explicitly why even 
though 0'^ (or tt) enters quadratically in the action, 
its equation of motion is not independent of the other 
Stiickelberg fields, so not all of them are dynamical, and 
the BD ghost is hence absent. Now to go further we show 
how the same result remains true when all Stiickelberg 
fields arc excited, and consider an arbitrary potential, 
[11] 



,(2) 



,(3) 



Cm = ^^er(^) + "3/:^::(/C) + a.C'Zi^) (20) 



(4), 



with /C(^ = (5^ - ^/di^i)°^d^AFn^h and 
4^2 = [/C]3-3[/C][/C2]+2[/C3], 

'(4) 



(21) 



4cr = [/C]4-6[/Cl[/C]" + 8[/C-^][/C]+3[/C^]"-6[/C4]. 

We then define the associated energy-momentum tensor 
as 



X„ 



f^fjai/ '^ui' 



(22) 



+ (1 + 3a3) (iCl, - ICIC^u + I ([/C]2 - [/C^]) g,,^ 
+ (aa + 4^4) (iCl, - ICJC% + i/C^, {[ICf - [JC^]) 



-([/C]3-3[/C][/C2]+2[/C3]).9^, 



such that the equations of motion are easily expressed as 
£a - a^A^ = . (23) 

The kinetic matrix is then 

9A" 5A0 



<^ab ~~ 



(24) 



To prove the absence of ghost, it is sufficient to show that 
det A = 0. 



Fourth-order expansion 

For simplicity, we write in this section the Stiickelberg 
fields as (p"^ = x'^ + V^ and address the question order 
by order in powers of V"'. We will stop the analysis 
at quartic order in V^ (in the Lagrangian, eubie order 
in the equations of motion) as it is at that order that 
potential problematic terms have been found previously. 
Once again we work in the locally inertial frame (as we 
show later, including gravity bears no effect on the re- 
sult). Setting a^ — a^ ~ 0, the equations of motion are 
then 



is unclear whether such a background is stable and could 
be generated naturally in a physical system, as is ex- 
plained in [18]. In particular we know from the ADM 
formalism, that the shift cannot be excited without si- 
multaneously exciting the 3d metric, [11, 12]. However, 
for the purpose of this work we follow here the same pro- 
cedure as [17] and expand around this given background. 
We therefore focus on the quadratic fluctuations of the 
Stiickelberg field around this metric configuration, and 
to start with, we only keep track of the time derivatives 
(spatial derivatives are introduced later). 



£ -^^' ^7^(1) + J^(2) , x-(3) 



0. 



(25) 



Primary Constraint 



where J^i" is n*^ order in V^: 

J"(^' - J" 



(26) 






+ ^:F^^d^v^di3V" + ii^, 



^{dcV^f 



4 8 



1 

12 



d'^V^do^V-ydpVc 



(27) 



with E^ix = d^Vi, + di,V^. Here again, while neither 
J^g ' ' depend on (f)'^ , the next order correction J^^J has 
an explicit dependence on (jf' . However one can ex- 
plicitly compute the determinant of the kinetic matrix 
which turns out to vanish to that order in the expan- 
sion, 1^1 = 0{{dVY), so once again, only three of the 
four Stiickelberg scalars are dynamical. In this case, the 
associated constraint is expressed as 

C = £o + (\:F^^ - \^l^ + \^o,Va^^^^"\ £, - , (28) 



The potential is then of the form, (to all orders in t) 



i^ri 



2VT^ 



Vf.uV^'V-' + 






(29) 



\^/^ + £il/V° + 2ldjV'V^ + \il{V'^f + 



and it is clear from the last term (\{y^Y that V^ seems 
to acquire a kinetic term. However the Hessian matrix 
associated with the kinetic terms of the Stiickelberg in 
(29) is clearly degenerate. 






det yt = , (30) 



which can be easily seen by setting W^ ~ V'^+l^V'^, such 
that the Lagrangian only depends on W^, 



r — -- 



1 



271^7 



^^2 



{W 



{irW) 



^-^l 



(31) 



and it is clear that only three of the four Stiickelberg 
fields are dynamical. We emphasize however that this 
change of variable is not required to count the number of 
dynamical degrees of freedom or to see the constraint. 



C = PQ-tP, = ^, 



(32) 



and involves no double time derivatives. 



where P^ is the conjugate momentum associated with 
y^, and we neglect spatial derivatives. 



Non-trivial background 

It is clear from the previous expansion, that V^ seem to 
acquire a kinetic term at quartic order in the expansion, 
or in other words, when expanding around a non-trivial 
background. To make this more explicit, and to connect 
with [17], we give up the LIF in what follows and con- 
sider a specific example where the metric has a non-trivial 
background. Namely, we consider the following configu- 
ration, ds2 = -df2 -H 5ij- (dx* -h 2Pd£) (dx^ ^ 2£Jdi). It 



Secondary Constraint 

So far, we have neglected any space-like derivatives as 
the previous analysis was sufficient to prove the existence 
of a primary constraint, and hence already remove half a 
degree of freedom, already implying the absence of ghost. 
However to study the algebra of the constraint system, 
and the existence of a secondary constraint, it is impor- 
tant to re-include space-like derivatives. When doing so. 



one can easily check that the canonical momenta associ- 
ated with the fields y are then 

Po = i'J'^0-ilV" + 2tf^^V,-£l^^V' (33) 

p, = j-,0 - e^v" - 2eAv'' + \td,Vk (34) 



1 



fcT> 



where for simplicity we stop the expansion at quadratic 
order in t. It is then straightforward to check that the 
conjugate momenta are not independent but are related 
in the following way: 

C = Po- PP., - £| d^V + ted^Vj = 0{t) . (35) 

This relation, can still be regarded as a primary con- 
straint as it appears at the stage of determining the mo- 
menta. Furthermore, requiring the conservation of the 
constraint (35), leads to the secondary constraint C2, 



C2 



D 



{C,U} 
d,P' + 



= 



(36) 



id,P' 



To summarize, the Hamiltonian is accompanied by the 
two constraints, (35) and (36). These two can be used 
to eliminate one of the coordinate-momentum pair (for 
instance (F^, P3), although any other pair is equivalently 
acceptable). Hence, the system propagates three phys- 
ical degrees of freedom. In particular, from the pri- 
mary constraint (35) one can express P3 in terms of 
Po, Pi, P2,V^ ,V^ ,V^ , and from the secondary constraint 
(36) one can express V^ in terms of Pi,P2,V^,V^,V^. 
Hence, the physical Hamiltonian depends only on the un- 
constrained variables: 



fl-phys fl- 



phys 



iPo,Pi,P2,V'>,V\V^). 



(37) 



OUTLOOK 



The theory of massive gravity as proposed in [11] has 
been shown to be free of ghost both in the decoupling 
limit [10] and to fourth order in perturbations in the 
ADM formalism. Recently, Hassan and Rosen have ad- 
vanced the analysis further and shown the absence of 
ghosts to all orders in the ADM language, [12]. While 
this result is in complete consistency with what was found 
in [10, 11] recent works seem to suggest the existence of a 
ghost cither at quartic [17], or even worse at cubic order, 
[19] in the Stiickelberg language, [20]. If the conclusions 
of these works were correct one should seriously investi- 
gate why the ADM and Stiickelberg language lead to so 
different answers. In this work we have filled the missing 
gap between both languages and shown the existence of 



a constraint among the four Stiickelberg fields which is 
responsible for the absence of ghosts. We perform this 
analysis to all orders in specific toy-models, as well as to 
quartic order in the full-fledged theory. In all these cases 
we show that the Stiickelberg analysis is in complete con- 
sistency with the ADM language [11, 12] and proves the 
absence of the BD ghost in this family of massive gravity 
theories. 
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To be precise, Ref. [19] uses a similar but not identi- 
cal methodology to diagnose the presence of a ghost by 
looking for higher derivative interactions of the helicity 
zero mode. Although the present article is concerned with 
Stiickelberg decomposition, it is clear that for similar rea- 
sons the analysis of [19] is incomplete, since within the 
helicity formalism, extra constraints will arise (as is in- 
evitable for consistency with the ADM formalism), and 
these constraints will be seen to remove the apparent 
ghost. 



